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00 ■ Abstract 

Stochastic differential equations for processes with values in Hilbert spaces are now largely 
^v^j | used in the quantum theory of open systems. In this work we present a class of such equations 

■ and discuss their main properties; moreover, we explain how they are derived from purely 
quantum mechanical models, where the dynamics is represented by a unitary evolution in 
a Hilbert space, and how they are related to the theory of continual measurements. An 
essential tool is an isomorphism between the bosonic Fock space and the Wiener space, 
which allows to connect certain quantum objects with probabilistic ones. 



o 



OS 



1 Classical stochastic differential equations 

Recently there was an increasing use of stochastic differential equations (SDEs) in quantum 
optics and more generally in the theory of quantum open systems |l ]— 1|. In particular SDEs 
are used for theoretical purposes (2|, [3|, ||] in the theory of quantum continuous measurements 
H f?J, when a system is continually monitored in time, and for numerical simulations of quantum 
master equations [Q], which are mathematical models for the dynamics of quantum open systems. 
In this work we want to describe a restricted class of such equations, to discuss their main 
properties and to show how they can be derived from fully quantum mechanical models ||] . 

We consider only SDEs of diffusive type, when the noises are Wiener processes; for the theory 
of Ito's stochastic integrals in Hilbert spaces we refer to |[(]]. Let [Q, (Ft), J 7 , P) be a stochastic 
basis satisfying the usual hypotheses; in particular, (fi,./ 7 , P) is a probability space, which we 
assume to be standard, (^)t>o is a filtration, i.e. Tt is a sub-a-algebra of T and T s C T t 
for s < t; we also assume T = VtX)^ 7 *- Let Wfc(-), k = 0, 1, . . . , be a sequence of continuous 
versions of adapted, standard, independent, Wiener processes with increments independent from 
the past. 

Let Ti be a separable complex Hilbert space and let us consider the following linear SDE 
with "multiplicative noise" for an "H-valued process tpt- 

oo 

d^ = ^RjWikdWM-iKWikM 

3=0 ^ L - L) 

ipo = ip 
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where the coefficients Rj(t), K(t) are operators in TL satisfying some assumptions to be discussed 
below and the initial condition ip is an W-valued .Fo-measurable random variable. Let us state 
now some different sets of assumptions about the coefficients. 



Assumption 1.1 The coefficients are time dependent, bounded operators on TL satisfying: 

(i) ^2^L Q R*jif)Rj{t) is strongly convergent in C(Tt) for all t; 

(ii) the functions t i— > R(t), t i— > -ftT(i) are strongly measurable, where R(t) is the operator from 
TL toTL®l 2 defined by (x®c\R(t)y) = Y%L Cj (x \ Rj (t)y), Vx,y G TL, Vc G Z 2 ; 

2 



(iii) VT 6 



ess sup 

t<T 



j=0 



+ \\K 



< oo . 



Assumption 1.2 TTie coefficients are time independent, (in general) unbounded operators, de- 
fined on a dense domain T> C TL; we assume that 

oo 

\\Rp\\ 2 n < oo, VxeP*, 

j=0 

where T>* is a core for the adjoint operators K* and R*. Moreover we assume the following 
dissipativity condition to hold: 

oo 

Yl W R 3 x fn + 2Im ( x \ Kx ) n - c W x fn ' Vx G £>. 
j=0 



Assumption 1.3 The coefficients are time independent, (in general) unbounded operators, de- 
fined on a dense domain V C TL, satisfying 



•X- 



< oo , Vx G £>; 

j=0 

moreover, —\K is the generator of an analytic operator semigroup T(t). 



Let us remark that, in the time independent case Assumption 1.1 is stronger than Assumption 



1.2 and indeed Theorem |1 . 1| will show that, under Assumption 1.1, we can obtain stronger results 
than under Assumption |1.2j . We note also that, for a linear equation, the so called "Lipschitz 
conditions" pi are equivalent to the boundness of the coefficients. 



Definition 1.1 We consider three different kinds of solution: 

(i) a (topological) strong solution of eq. (|1.1[) is an "H-valued process such that Vt G M+ the 
following equation holds almost surely (a.s.) 

^ = V + y"/ R J (s)iP s dW j (s)-i K(s)iP s ds; (1.2) 

j=0 Jo Jo 
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(ii) a (topological) weak solution of eq. ( |1.1| ) is an 7Y-valued process such that Vi E R+ and 
Vx G T>* the following equation holds a.s. 

(x|&>„ = W>„+Z;/ (R* j ( s )x\^s) H dW j (s)-i (K*(s) X \iP s ) n ds; (1.3) 

•_ Q J -/ 

(hi) a (topological) mi/rf solution of eq. ([hi]) with Assumption L3 is an W-valued process such 
that \/t E R+ the following equation holds a.s. 

00 ,. t 

ih = T(t)iP + V / T(t- s)Rjip 8 dWjis) , (1.4) 
j=0 Jo 

where T(s) is the analytic semigroup operator generated by —\K. 

For more details about these definitions see |1C]. It is clear that a strong solution is also a 
weak solution and it is possible to show that a weak solution is also a mild solution (see [p7cf| , 
Propositions 6.2, 6.3, 6.4 and Theorem 6.5). In the following theorem we state some results 



about the existence, the uniqueness and the Markov properties of the solution of eq. (1.1) under 
different hypotheses. 

Theorem 1.1 (i) Under Assumption eq. (|1.1| ) admits a unique strong solution; (ii) under 
Assumption eq. ( |1.1[ ) admits a unique weak solution; (iii) under Assumption eq. (1.1) 



admits a mild solution and this is also a weak solution. All these solutions are strong Markov 
processes. 



Proof. The statement (i) is Theorem 1.1 of [11|; for the proof of (ii) see Theorem 1 of ]12| and for 



the proof of (iii) see Theorems 6.19 and 6.22 of [10|. For the Markov property see [12| Theorem 



1, |H| Theorems 9.8 and 9.14. □ 

We are not interested in eq. ( |1.1[ ) in general, but only when HV'tll^ is a martingale and can be 
interpreted as a probability density with respect to P. Indeed, if we have a positive martingale 
HV'tll^ with E p [H^tll^] = 1 (Ep denotes the expectation with respect to P), then 

P{A) :=E P [l A Ut\\ 2 M ] , VAG^,ViGM+, (1.5) 

defines a new probability law P on (Q,^). To obtain this we need a further assumption. 

Assumption 1.4 For all t G M + , the coefficients satisfy the following operator equality 

00 

j=0 

(where the sum is convergent in the strong operator topology) and the initial condition is a square 
integrable, normalized random variable, i.e. E p [HV'II^] = 1- 

Theorem 1.2 Under Assumptions [7j| and the square norm WiptW^of the solution is a (pos- 
itive) martingale and E p [llV'tll?,] = 1- Moreover, under the law P defined by eq. (|i~E]), the 
processes 

W k (t):=W k (t)-2 [ Rem k {s)ds, k = 0,1,... , (1.6) 
J 
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where 



m fcW := iiTli? ' ( L7 ) 



IIWI„ 

are independent, adapted, standard Wiener processes. 



For the proof of this theorem see [11], Theorem 1.2 and Proposition 1.1. The unbounded case 



is more difficult and it is treated in [12] under an additional condition called of hyperdissipativity. 



Let us introduce now the process 

i>t ■= n /n exp f — i / (Imm k (s))(dW k (s) -Rem k (s)ds) S )'ipt- (1.8) 
\m\\n V f^Jo J 

Note that = lj indeed, exp(- • • ) is a stochastic phase with no particular meaning, which 

we introduce only in order to simplify the form of the stochastic differential of ipt- Note also 
that fh k (t) = (ipt Rk(t)ipt)- 



Theorem 1.3 Under Assumptions |7T1] and \1.4 , in the stochastic basis ($7, (Ft)i P), ipt sat- 
isfies the non linear SDE 



oo 

d^ t = -ik (t, dt + Rk (t, dw k (t), (u 



k=0 



where, V/ G H, 



K(t, f) := \ m) + K*(t))f - i ± (^|^ W) ~ *(*)) f + 



2 h \ k{) — ) \ k{ ) — mr ) f • ( 0) 

and 



(f\L k (t)f) , 

/?,(/■/):- ( ) /. (.1.11) 



Moreover ipt is the unique strong solution of eq. ( |1.9| ) and it is a strong Markov process 



Proof. The first part of this theorem is essentially Theorem 1.3 of 11]. Moreover it is easy 
to show that the coefficients K and R k satisfy the local Lipschitz conditions (see p. 198); 
hence ipt is the unique strong solution of eq. (11.91). Again the Markov property follows from [fix]] 



Theorems 9.8 and 9.14. □ 

Let us stress that the Markov property allows to associate to the process ipt a transition 
function satisfying a Chapman-Kolmogorov equation (see [10] Section 9.2.1 and, in particular, 



Corollary 9.9); similar considerations apply to the solution of the linear SDE. 

Examples of equations of the types (1.1) and ( |1.9| ) were introduced in the theory of quantum 



continual measurement in Ref. 
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In the statistical formulation of quantum mechanics the states of a quantum system are 
represented by statistical operators. Let T(H) be the trace class on H and S(Tt) denote the set 
of statistical operators, i.e. 

T(H) := {g G £(H) : Tr n {^} < 00} , 
S(H) := {g G T(H) :q* = Q,Q>0, Tr n {g} = 1} . 
Equipped with the norm \\g\\i := Tr n { \J g* g } , T(7i) becomes a Banach space whose dual is 



Theorem 1.4 Let tpt be the solution of eq. ([O]) with Assumptions LI and 1.4 - Then, the 
formula 



Ti n {g t a} = E p {ip t \aip 



Va G C(H), 



(1.12) 



defines a family {gt, t > 0} of statistical operators on TC with go = g, where Tr^ {ga} :- 
■p\aip) n , Va G CiTL); eq. ( |1 . 12 ) is equivalently written as 



E 



Tr H {ot a} = E 



Va G £(H). 



Moreover, gt is the unique solution of the integral equation 

g t = g+ / £(s)*[g s ] ds, 
Jo 

where >C(s)* is the preadjoint operator of £{s) G C(£(7i)) defined by 

1 oo 

C(t)[a] := -~[K(t) + K*(t),a} + ~J2{Ri(t)[a,Ri(-t)) + [R* (t) , a] R^t)) 



(1.13) 



(1.14) 



i=0 



Va G 

(1.15) 



where [a, b] := ab — ba. 



For the proof of this theorem see [11 1 Theorem 1.2 and Section 3. 



In the theory of open quantum systems eq. (1.14) is usually written as dgt/dt = C(t)*[@t] (if 
C(t) is sufficiently smooth in time), it is called a master equation and describes the evolution of 
the statistical operator representing the state of the system. 

Formula ( 1.13| ) connects the solution ipt °f the SDE ( |1.9|) with the solution gt of the master 
equation (1.14); this fact is at the basis of stochastic methods for numerical simulations of 
solutions of master equations [Q]. 

For what concerns the theory of continual measurements, we shall see in Section 3 that 
the processes Wk(t), or functionals of them, represent the output of the measurement and 
that the physical law of this output is given by the probability P defined in eq. (|1.5j ). In 
particular, formulae for all the moments (under P) of the processes Wf.{t) can be given. For 
instance, from Theorems 1.2 and 1.4 we get E_ [W^(t)] = 2 J * ReE. [m&(s)] ds and E_ [m^(s)] 



Tx H {g s Rk(s)}; therefore, the mean value of the output Wfc(i) is connected to the solution of 
the master equation by 



% [W k (t)\ = fTr H {g s (R k (s) + R* k (s))}d S . 
Jo 



(1.16) 
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The computation of the moments of the second order is a little bit more involved; one needs 
to go through the characteristic functional |7|, 13] or to exploit the Markov properties of the 
solution of eq. (|1.9|) . Let U(t,s) be the evolution operator (from time s to t) associated with 
eq. ( PH) , i.e. g t =U(t,0)[g]; then, the result is (cf @ eq. (5.20)) 



E p [W i {t 1 )W j {t 2 )]=t 1 M 2 + 

+ j dsi j ds 2 Ti n |^i?j(si) +Ri(s 1 fjU(si,s 2 ) [Rj(s 2 )g S2 + Q S2 R*j(s 2 )]} 

+ J ds 2 J^ dailt w {(i^(a 2 )+iJ5(fl2))W(fl 2 ,«i)[i^(ai)e. 1 + e. 1 ^(«i)]}. (1-17) 



+ 



Also the random vectors Vt have a physical interpretation 0, ||, [|, Indeed, it can be 



shown that the random orthogonal projection 



represents the state at time t one has 



to attribute to the quantum system if the trajectory of all the processes Wj is known up to 



time t; in this sense 



Qt 



E, 



represents the a posteriori state of the system at time t. Then, 
is the state of the system when no information on the outputs is taken into 



account; therefore, Qt is called the a priori state at time t. 



Note that no physical quantity, such as Qt, 



and P, depends on the arbitrary phase 



introduced in eq. (l.< 



2 The Fock and Wiener spaces 

In order to be able to show how the stochastic equations of the previous section are related to 
the quantum dynamics, we need to recall some facts about Fock and Wiener spaces. 

Let Z be a separable complex Hilbert space and let us denote by V the symmetric Fock space 
over X := Z ® L 2 (R) ~ L 2 (R; Z), which is defined by 

oo e 

r:=c ©E (<^ n )s> 

n=l 

where {X ® n ) s denotes the symmetric part of the tensor product X® n and is called the n-particle 

space (0 Sects. 17 and 19). The vectors E{f) := (l, /, / ® /,..., -±= /®» ...),/€*, 

are called exponential vectors (Jj] p. 124); let us denote by £ the linear span of the exponential 
vectors. 

Theorem 2.1 The set of all the exponential vectors is linearly independent and total in T; 
moreover, for f,g£Xwe have (E(g)\E(f)) r = exp(g\f) x . Let us fix g £ X and V £ U(X) 
(unitary operators on X); there exists a unique unitary operator W(g;V) (Weyl operator) such 
that 

W(g;V)E(f):=exp{~\\gf x - (g\Vf) x \ E(Vf + g), V/ € X. (2.2) 

These operators satisfy the multiplication rules fc/j £ X , Vi G 14{X)\ 

W(g 2 ;V 2 )W(gr,V 1 ) = exp[ikrL(V 29l \g 2 ) x }W(V 2 g 1 +g 2 -, V 2 V X ) . (2.3) 
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Total means that £ is dense in T; for the proof see [14] p. 124, Proposition 19.4 p. 126, pp. 134- 
135. Note that = exp||/||^,; by normalizing the exponential vectors one obtains the 

coherent vectors e(f) := exp ] E(f). The vector e(0) = E(0) = (1,0,0,...) is called 

the vacuum. 

In the following we shall need particular Weyl operators. Let V G IA{X) be such that 
(Vf)(t) = V(t)f(t), V/ G X ~ L 2 (R;Z); then V(t) G U(Z). Let g G ^ oc ( R ; 2 )> 
Is I|3( T )II| dr < +oo, Vs,i G R, s < i. Then, we set 

W t (g; V) := W ( X[0) t]<7; X[o,t]( V " + *) ■ ( 2A ) 

The Fock space has a continuous tensor product structure, i.e. Vs,i G R, s < t, T = VL^ 
r* (g) r^"°°, where is the symmetric Fock space over L 2 ((a, 6); i?) , — oo < a < b < +oo. 
Correspondingly one has E(f) = E (x{-oo,s)f) ® E [x( s ,t)f) ® E (x(t,+oo)/); w e also denote by 

C T b a the linear span of the exponential vectors of the type E (x(a,b)') ( Fl Proposition 19.6 
p. 127 and pp. 179-180). Note that W t (g; V) leaves invariant 1%. By setting S{t)E(f) := E(f t ), 
ft{s) := f(s — t), one defines a strongly continuous unitary group of shift operators on T; 
moreover, one has S(t)T^ = T^ t . 

Let {ej, j = 0, 1, . . . } be a c.o.n.s. in Z, g G X and let us define on £ the operators 



a(g)E(f):=(g\f) x E(f), V/ G X, A^t) := a (x[o,t]U e j) > ( 2 ' 5 ) 

d 
de 



a\g)E(f):=^E(f + eg)\ £=o , V/ G := a* (x^O)^) J (2-6) 



let us note that one has {a) {g)E{fi)\E(f2)) T = {E(fi)\a(g)E{f2)) r . The operator a){g) turns 
out to be linear in its argument and a(g) antilinear. The domains of these operators can be 
suitably extended; there exists a core V for all these operators where all their products are 
well defined and such that T> D £ , T> D (X® n ) s for every n. Then, Aj{t) and a{g) send the 
n-particle space into the (n — l)-one and Aj(t) and o){g) send the n-particle space into the 
(n + l)-one; by this they are called annihilation and creation operators, respectively. On T> 
the annihilation and creation operators satisfy the canonical commutation rules [a(g), a(f)] = 0, 
[a' (5), = 0, \a(g), «'(/)] = (g\f) x -> which characterize Bose fields (H| pp. 136-146 and 

Example 24.1 p. 181). For the operators Aj(t) and Aj(t) these commutations relations become 
[Aj(t),Ai(8)] = 0, [A](t),A\(s)} = 0, [Aj(t),Al(a)] = %min(i,s); in particular 

[A j (t) + A](t),A i (s) + Al(s)] = 0, Wi,j,t,s. (2.7) 

A connection between creation and annihilation operators and Weyl operators is given by 
the following theorem (0 pp. 136, 139, 143-144). 

Theorem 2.2 For every k G X the operator a(k) + a\k) is essentially selfadjoint on £. 
Moreover, if we denote again by a(k) + a){k) its selfadjoint extension, we have W(ifc; 11) = 
exp{i[a(k) +a)(k)}}. 

^From now on, when v G Z, we shall set Vj := {ej\v)z- Moreover, we shall confuse essentially 
selfadjoint operators with their selfadjoint extensions. 

Let us set K := {k G Lf oc (R + ;Z) : kj{t) G R, Vj,*}. From eq. (§]| we have that 



[W tl (ifc i; 1) , W t2 (ik 2 ; t)] = , Vti, t 2 G R+ , Vfc, fc 2 G ft, (2. 



7 



or, which is the same, from the canonical commutation rules we have 



(X[o,ti] fc i) + flt (X[o,ti]&i) > a (X[0,t2} k z) + flt (X[0,t 3 ]*2) = • ( 2 - 9 ) 

Then, we set ,M(t) := {Wt(ifc; 11), fc G 72}" (double commutant in C(T)); M{t) is the von 
Neumann sub-algebra of C(T) generated by the Weyl operators {Wt(iA; 1), k G 72.} and it turns 
out to be a commutative algebra. Moreover, by the spectral theorem, the previous one and the 
definitions (|2.5| ) and (|2.6| ) we have also that every operator in M(t) can be seen as a bounded 
function of the commuting selfadjoint operators {Aj(s) + A^(s) , j = 0, 1, ... , s£ [0,t]}. 

Let us consider now infinitely many independent standard Wiener processes Wj(t) as in 
Section 1, but let us choose a concrete probability space: O is the set of all the functions 
uj '. t i — > uj(t) = (u>o(t) , uJi(t) , . . . ) with u!j(0) = and t \— * ujj(t) real and continuous, T is the 
o"-algebra generated by all canonical projections, P is the Wiener measure on (fi, T), the Wiener 
processes are given by the canonical projections, i.e. Wj(t)(u) = uij(t), Tt is the sub-cr-algebra 
of T generated by the family of random variables \ Wj(s), j = 0, 1, . . . , s G [0, i]}. Let us call 
L 2 := L 2 (Q,T, P) the Wiener space; we recall that for a,b G L 2 one has 



(a\b) Llr = / a(u;) 6(a;)P(cL;) = E p [a 6] . (2.10) 
'fi 



Theorem 2.3 There exists a unique unitary isomorphism I : Tq — ► L 2 satisfying 



f 00 r 

lE(f) = E w (f) :=exp ^ / 

Moreover, we have 

XWtiik-l)!- 1 

and, Vj = 0, 1, . . . , Vi > 0, 



f j {t)dW j {t)--{f j {t)fdt 



exp{if^\( S )dI^( S )}, 



V/GL 2 (M+;Z). (2.11) 



VA G 72, Vt > 0, 



l(A j (t) + AUtj)l- 1 = W j (t), ZMit)!- 1 =L°°(n,F t ,P); 



here Wj(t) is a multiplication operator in L 2 . 



(2.12) 



(2.13) 



Proof. The isomorphism X is constructed in [14] Example 19.9 pp. 130-131. 
By eqs. @, Q, ( gH| ), we have VA; G Z) and V/ G L 2 (M+; 2) 

JW t (i*;l) T- x E w {f) = expjif^ dWj(s) + (2f J (s)+ik j (s))lmk ] (s)ds\ }e w (J). 



(2.14) 

Then, Theorems [DJ U and eq. ( gjjj ) give eq. ( glp and the first of eqs. <^Wj ). By eq. (pl^ ) 
and the definition of «M(f) we have also the second of eqs. (2.13). □ 

It is possible to extend the isomorphism I to the whole T = T ^ (g) To by mapping it into 
l? w ® I? w . Let us stress that X is the unitary transformation which simultaneously diagonalizes 

all the commuting selfadjoint oper ators {Aj(t) + A\(t), j = 0, 1, . . . , t G 



8 



3 Quantum stochastic differential equations 



Let 7i be a separable complex Hilbert space as in Section 1. Hudson and Parthasarathy [ |l5| ] 
developed a quantum stochastic calculus which gives meaning to integrals with respect to dAj(t), 
dAj(t), dt- we give only a rough idea of the definition, just what we shall need in the following. 
Let G(t) be an operator with domain including TL(&6 (the linear span of the vectors of the type 
(p®E{f)) and such that G(t)H®£t 00 C H^TL^. Then, the integral f* G(s)dA,-(s) is defined 
as a suitable limit (if it exists) of a sum of terms of the type Gifj) (Aj(ti+i) — Aj(ti)); this kind 
of definition is analogous to the definition of classical Ito's stochastic integral (|l4[] Sects. 24 and 
25). For instance, one has 

00 r+oo 00 />+oo 

a(k) = J2 kJiFjdAjit), a \k)=^2 %(*) d ^(*)> k £ L 2 (R + ; Z). (3.1) 

•_ n JO -_ n JO 



3=0 



3=0 



Let us introduce now the operators H, Lj G £{7~t), such that H = H* and X]j=o LjLj 
strongly convergent in C(TC); let us set K := H — ^Yl'jLo L*Lj. We consider the quantum 
stochastic Schrddinger equation 



dU(t,s) 



oo 

E [ L id^(t) - L*dAj{t) - iKdt \u{t, s) 



U(s,s) = l; 



(3.2) 
(3.3) 



here s < t, U(t, s) is an operator on W§T and Lj is identified with Lj <g> 1, Aj(t) with 1 £g> 
and so on. 



Theorem 3.1 (Hudson, Parthasarathy, Frigerio) Equations ( \3.3J have a unique solution 

U(t,s), t > s; we have also U(t,s) G U(H g) T*) C U{H®T) and U(t,r)U(r,s) = U(t,s) for 
s < r < t. If we fix f G X and set Va G £(W) 



C t [a] := i 



X -i . ' 

H + i E (OT L i " £ > a + o E ( t L 3 » °\ ^ + L 3 [«« L 3'] ) » 



3=0 



(3.4) 



3=0 



then we have, V£ £ W, Va £ Vt, s £ 1, t > s, 

<*7(i, s)£ ® e(/)|a E7(t, s)£ ® e(/)) W8r = <£K>„ + 



(U(t, s)S e(f)\£r[a] U(r, s)£ ® e(/))„_ r dr. (3.5) 



Moreover, the quantity S*(t + s)U(t + s, s)S(s) does not depend on s and {U(t), t G K}, where 
U(t) := S*(t)U(t,0) fort>0 and U(t) := U*(-t,0)S(-t) for t < 0, is a strongly continuous 
one-parameter group of unitary operators on H (g) T. 

For the proof of the first part of the theorem see H] Theorem 27.8 p. 228, Proposition 



26.7 p. 216, Corollaries 27.9 and 27.10 p. 230; for the statements about U(t) see [16]. Let us 
stress that from the strong continuity we have that U(t) is strongly differentiable on some dense 
domain and the same is true for S(t). But the two infinitesimal generators are unbounded and 
mutually related in such a way that U (t, s) = S(t)U (t — s)S*(s) is not differentiable in the usual 
sense; however, although the usual derivative of U(t,s) does not exist, the quantum stochastic 
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calculus gives meaning to its stochastic differential (p.2|). Another interesting point is that, from 
the connection between U(t) and the master equation (|3,5| ), it is possible to deduce that the 
Stone generator of U(t) is a selfadjoint operator, necessarily unbounded from above and from 
below. 

In the standard formulation of quantum mechanics it is assumed that the evolution of an 
isolated system is represented by a strongly continuous group of unitary operators. So, by 
the previous theorem, we can think we have a quantum system A, described in the Hilbert 
space TC, interacting with a second system B, a bosonic field described in T. The dynamics of 
A + B is given by U(t); if we interpret S(t) as the free dynamics of the field, then U(t,s) is 
the evolution operator of the compound system in the interaction picture with respect to the 



free-field dynamics. It is possible to show [O, 13| that evolutions such as U(t,s) are a sensible 



approximation to the dynamics of a photo-emissive source A, such as an atom, an optical system 
in a cavity, . . . , interacting with the electromagnetic field B. 

Let us fix now s = as initial time and £ e(f) as initial state, £ G TC, = 1, 

/ G Lf oc (K; Z); we have not defined coherent vectors for such an /, but we shall see in the 
following that in our applications the queue of / does not matter and that it is often useful to 
consider periodic functions. Indeed the choice f(t) ~ exp(— iuiot) could represent a stimulating 
monochromatic laser. Moreover, we set Ut := U(t, 0), so that Ut£,<3 e(/) is the state of our 
compound system at time t. 

Up to here we have studied the evolution of the system. Let us consider now a measurement 
process in which some observables of the system are monitored with continuity in time; the 
mathematical model of the detection scheme will be given again in terms of operators on the 
Fock space and quantum stochastic integrals. Quantum stochastic calculus has been introduced 
into the theory of continual measurements in [O] . 

Let us recall that, in the usual formulation of quantum mechanics, observables are represented 
by selfadjoint operators. Let us consider a quantum system represented in a Hilbert space f) 
and let ipt be the state of the system at time t (ipt £ f), HV'tllf) = !)• Let X\, X n be 
commuting selfadjoint operators with joint spectral projections 11^(5), B Borel subset of W 1 . 
Then, (^|IL»(B)^){) is interpreted as the probability that the measurement of the "observable" 
X gives a result in B at time t; indeed, (tpt I IL^ (-)ipt )f> is a probability measure on R n . Obviously, 
the whole probability measure is known if we know its Fourier transform (the characteristic 
function), which turns out to be given by (^|exp{ifc • X}^)^, k G M. n . We can also say 
that we know completely the probability law of X if we know all the quantities {ijjt\a ipt) t) with 
a belonging to the commutative von Neumann sub-algebra of £(f)) generated by the set of 
operators {expfifc • X], k G M n }. This can be generalized even to families of infinitely many 
commuting selfadjoint operators. 

Now, let us consider again our system A + B and the family of commuting selfadjoint 
operators 

z iV)-=Y, V ij (s)dA i (s) + V ij (s)dAl(s) , t>0, j = 0, 1, ... ; (3.6) 
i=o J ° L J 



here V, V(t) are unitary operators as in eq. (2.4) and Vij(t) := {ej\V(t)ei} z . We consider {Zj(t)} 



as a set of continuously monitored observables. Indeed, in the case of electromagnetic field 
emitted by some source, the monitoring of observables of the type ( [3.6D can be concretely realized 
by a measurement scheme called balanced heterodyne detection [[b|. Let us denote by M.z(t) the 
von Neumann algebra generated by the Weyl operators jexp i^jlo Jo %( s ) dZj(s) , k G 72.|. 
By the previous discussion, the law of the measurement of Z(s) for s from to t is known if we 
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know the quantities (Ut£ <8> e(f)\YUt£, (g) e(/)) H0rQ i VY G M.z(t). The next step is to transform 
these quantities and to link them to the stochastic processes introduced in Section 1. 

^From now on V is the unitary operator appearing in the observables (|3.6|), / is the function 
characterizing the initial state of B, £ is the initial state of A (£ G Tt, \\^\\ n = 1) and Hq = 
Hq G £(T~L) is a selfadjoint operator which can be chosen arbitrarily and is introduced for future 
convenience. 



Theorem 3.2 For every t G M+ and Y G M.z{t), we set 

U t := e X p[i# i] W f *(/; V)U t W t (f; V) , 



(3.7) 



tk :=lU t £®e(0), 



(3.8) 



Y:=W*(f;V)YW t (f;V). 
Then, we have U = 1, Y G .M(i), TYX" 1 G L 00 ^,^, P), 
m ® e(/)|me ® e(/)> W9ro = ® e(0)|y^ e(O)) H0rQ 



(3.9) 



(3.10) 



d&t = lJ2 [^)d^(t) - i?*(t)d^(t)] - iK(t)dt \u t , (3.11 



(3.12) 



i=0 



(3.13) 



finally, ipt satisfies the SDE (1.1) wf/i initial condition Tpo(uj) = Vw G 12. Moreover, if we 
assume ess sup||/(i)|| z < +oo for all T G M+, i/ien i/te coefficients ( 3.12 ) and ( 3.13] ) satisfy 



t<T 



Assumptions 1.1 and 1.4 



Proof. Note that from eqs. fl2.2|) -(pO|) we have 

Wt(f;V)e(0)=e{ Xm f), W ((l - X[o,t]) /; t)W t {f;V)e{Q) = e{f) . (3.14) 

The operator Ut acts not trivially only oeH® Tq, while W ((l — X[o,t l) /> l) om y on F° <8> I\; 
therefore, these two operators commute. Then, by eqs. ( p.!4|) and ( |3.7|) one has 



ft £ ® e(f) = W ((1 - X [o,t]) /; 1) Wt(/; y)e- lH »* C/ t £ ® e(0) . 



(3.15) 



The first of eqs. ( |3.10| ) follows from eqs. ( |3.15| ), (|3.9| ) and the second one from eq. ( |3.8|) . The 
fact that Y belongs to M(t) follows from the definitions of M(t) and Mz{t) and from eq s. (|2.3| ) 
and (0); the fact that one has XYX~ l G L°°(r2,jr t ,P) follows from the second of eqs. fl2.13|) . 
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/From eqs. Q, (gj), (gj), one obtains how the operation W t *(/;V) • W t (/;V) 
transforms the operators Aj(t), A^(t); then, eqs. ( |3,11| )-( |3.13| ) follow from eqs. ( |3.7| ), (|3.2|). 

By eq. ( |2.5| ), the operators Aj(t) annihilate the vacuum. Moreover, by the factorization 
properties of T and the definition of quantum stochastic integrals, we have that dAj(t), acting 
on r* +dt , commutes with Ut, acting on H Tq, or, better, f Q ■ ■ ■ dAj(s)U s = L ■ ■ ■ U s dAj(s). 
Therefore, when eq. ( |3.11 ) is applied to £<8>e(0), the integral with respect to Aj gives a vanishing 



contribution and the integrand can be changed at will; so we can write IE] 



dU t £ ® e(0) = \ Rj(t) d A ](t) + dA j(t) 



3=0 



iK(t)dbjUt€®e(0) 



(3.16) 



By the first of eqs. ( |2.13j ) we have that ipt satisfies the SDE ( |1.1| ) with initial condition tpQ 
£ (g) £"^(0); but E w (0) = 1, so that we can write ipo(uj) = £, \/uj G $7. 



Finally, it is possible to prove that 



H o\\ + \ 




+2||/(t)|| 


'-J '■ > 


Assumptions |1.1| anc 


O. 





E^O R j (^j (*) - E^Lo L *j L 3 
1/2 



and II if 



< \\H- 



; then, one can check that our coefficients satisfy 

□ 



Let us comment the content of Theorem 3^. All the physical quantities (probabilities, char- 
acteristic functional, moments) are given by the "quantum expectations" (Uf£ <S> e(f)\YUt£ <8> 
e (/))w®r with Y G Mz(t). If we set 

Y := TY1- 1 = XWt(f\ V)YWt(f; V)!' 1 , (3.17) 

we have that Y is a random variable in (S7, Tt-, P) or (^l, J-t, P^j ; moreover, by the definitions of 
and P, we can write 

/. II 2 



Y^ t 



Y 



E. 



Y 



Finally, by eq. ( 3.10| ) we obtain 

(UtZ®e(f)\YU t ti®e(f)) n0ro 



E. 



Y 



(3.18) 



(3.19) 



where Y is given by eq. (p7|) , P by C5|) , is the solution of (Q with non random initial 
condition £, Rj(t) and if(i) are given by eqs. ( 3.12j ) and ( 3.13| ). Therefore, P is the physical law 
of the output, as stated at the end of Section 1. 

To discuss concrete physical applications of the previous formalism would be too long. How- 
ever, let us stress that the first non trivial example can be realized in the Hilbert space TL = C 2 . 
Indeed, in this space one can construct a model describing a two-level atom of resonance fre- 
quency u stimulated by a laser of frequency ojq and emitting fluorescence light; we are interested 
in the spectrum of the emitted light. Let us just give the list of the choices that character- 
ize the model: TL = C 2 , H = lo ^ with u > 0, Lj = aj ^ ^ with ay € C and 

< Ylj \ a j\ 2 < fj(t) = Aj e~ lw °* with ujq > 0, Xj £ C and Ao = 0; the arbitrary selfadjoint 

A N 



operator Hq is chosen to be Hq 



u 







(this eliminates some time dependence from the 



resulting equations). The fluorescence light is made to interfere with a strong laser of frequency 
v and the intensity of the resulting light is measured with a photodetector. It can be shown 
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that this scheme corresponds to a continual measurement of one of the components of Z (say 
Zq) given by eq. ( ^ ) with Vij(s) = e~ lus 5{j, v > 0. The intensity of the emitted light can be 
obtained from the second moments E p [W (ti)Wo(^)] of W (t) = V)Z (t)VV?(/; V)J _1 

; as a function of u, this intensity is the spectrum of the atom. Computations are given in pcf | 
where only quantum stochastic calculus is used (classical SDEs are not explicitly introduced); 
the result is a three-peaked spectrum, already known in the quantum-optical literature as the 
Mollow spectrum. 
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